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1.    BACKGROUND 

The  concentration  of  ozone  in  the  Toronto  area  exceeds 
Ontario's  regulatory  limit  more  often  than  any  other  pollutant. 
During  the  year  Ontario's  Air  Quality  Index  was  first  introduced, 
(June  1988  to  May  1989),  Shenfeld  and  Yap  (1989)  reported  that 
between  80  and  95%  of  the  hours  that  the  Air  Quality  Index  was  in 
the  Moderate  to  Poor  range  within  Metropolitan  Toronto,  ozone  was 
the  major  pollutant.   Furthermore,  during  the  summer  of  1988,  the 
ozone  hourly  averages  were  frequently  well  above  120  ppb, 
(AQI=50),  which  is  known  to  cause  an  impact  on  health  to  a  large 
percentage  of  the  population,  especially  those  who  may  be 
exercizing  outdoors.   These  high  ozone  developments  occurred 
during  periods  of  long-range  transport  (LRT)  of  ozone  and 
precursors  into  the  Province  from  the  United  States .   During  the 
Slimmer  of  1989,  the  occurrences  of  high  ozone  readings  were  much 
less  frequent,  as  the  LRT  contributions  were  not  as  prevalent, 
even  during  ozone  conducive  weather  such  as  warm,  sunny  and  light 
wind/calm  conditions . 

Kurtz  et  al.  (1989)  have  utilized  the  long-range  transport 
Acid  Deposition  and  Oxidant  Model  (A.D.O.M)  to  show  the  influence 
of  U.S.  emissions  of  nitrogen  oxides  and  reactive  hydrocarbons 
to  the  levels  of  ozone  in  Southern  Ontario.   Air  mass  back 
trajectory  studies,  most  recently  by  Yap  et  al.  (1988) 
established  that  U.S.  sources  of  these  pollutants  were 
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responsible  for  much  of  the  high  ozone  values  reported  in 
Southern  Ontario. 

Ozone  studies  by  Shenfeld  (1977),  Mukairanel  et  al.  (1982), 
and  more  recently  Yap  et  al.  (1988),  established  the 
meteorological  conditions  required  for  the  development  of  high 
ozone  levels  in  Southern  Ontario.   Studies  of  ozone  developments 
in  other  geographical  locations  by  Revlett  (1978),  Wolff  and  Lioy 
(1978),  Prior  et  al.  (1981),  Simpson  et  al.  (1983),  and  Robeson 
and  Steyn  (1990)  have  established  that  peak  hourly  average  levels 
of  ozone  can  be  correlated  to  specific  meteorological  parameters. 

2.    STUDY  OBJECTIVES 

The  main  objective  of  this  study  is  to  determine  the 
proportion  of  the  peak  ozone  levels  in  Toronto  resulting  from  the 
LRT  of  ozone  and  precursors  into  Toronto.   A  second  objective  is 
to  determine  whether  the  relationship  between  peak  ozone  levels 
and  meteorological  parameters  can  be  used  to  predict  the  peak 
ozone  levels  and  hence  the  Air  Quality  Index  for  the  afternoon 
and  early  evening.   The  third  objective  is  to  determine  the 
relationship  between  the  peak  hour  ozone  level  and  the  maximum 
8-hour  running  average  of  ozone  the  same  day  and  what  levels  of 
8-hour  averages  relate  to  peak  hour  levels  of  80  ppb  (Moderate 
air  quality)  and  120  ppb  (Poor  air  quality). 
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3.    STUDY  DESCRIPTION 

For  this  study,  we  have  selected  the  peak  hour  ozone  levels 
during  each  day  from  May  to  September  inclusive,  for  the  years 
1988  and  1989,  during  which  the  Pearson  International  Airport 
reported  temperatures  equal  to  or  exceeding  23° C  and  the  cloud 
cover  was  equal  or  less  than  four  tenths  opacity.   The  selected 
days  were  then  separated  into  days  that  LRT  occurred  from  the 
industrial  States  and  those  that  it  did  not.   This  classification 
was  determined  by  using  synoptic  weather  maps  and  the  Ontario 
Ministry  of  the  Environment  Air  Resources  Branch  archives  of  the 
48-hour  back  trajectories  of  air  into  Toronto. 

Three  Metropolitan  Toronto  ozone  recording  locations  were 
selected  for  study,  namely;  Downtown  Toronto  (Breadalbane) , 
Scarborough  (Lawrence  and  Kennedy)  and  North  York  (Yonge  and 
Finch) . 

For  LRT  and  non-LRT  days,  multiple  regression  equations  are 
used  to  relate  the  peak  ozone  levels  at  each  of  the  above 
locations  with  the  Pearson  International  Airport  maximum 
temperature,  average  wind  speed  and  average  cloud  cover  during 
the  hours  of  9AM  to  4PM.   Another  dependent  variable  used,  is  the 
peak  hour  ozone  level  measured  the  previous  day  at  the  rural 
station  outside  Simcoe,  Ontario.   The  difference  between  the  peak 
ozone  levels  deteirmined  from  equations  with  LRT  and  with  non-LRT 
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is  used  to  determine  the  contribution  of  ozone  to  the 
Metropolitan  Toronto  locations  due  to  LRT. 

4.    REGRESSION  ANALYSIS 

A  multiple  regression  analysis  was  completed  for  the  three 
sites  (Toronto,  Scarborough  and  North  York)  for  each  of  LRT  and 
non-LRT  days.   A  simple  regression  equation  was  derived  from  the 
ozone  data  set  of  the  Ontario  Ministry  of  the  Environment: 

Y  =  Bo  +  Bi  *Xi  +  B2  *X2  +  Bs  *X3  +  B4  *X4      (1) 

where:    Y   -  Maximum  Ozone  Level  (PPM) 
Xi  -  Maximum  Temperature  ( °  C ) 
X2  -  Average  Wind  Speed  (Km/hr) 
Xa  -  Average  Cloud  Opacity  (lO'^"»  ) 
X4  -  Previous  Day's  Maximum  Ozone  (PPM) 
Recorded  at  Simcoe,  Ontario 

The  results  of  the  regression  analysis  are  outlined  in  Table  1. 

The  partial  regression  coefficients  derived  from  equation 
(1)  are  not  independent  of  the  particular  metrics  used  for  the 
variables  and  are  therefore  not  directly  comparable.  Each  one 
gives  the  rate  of  change  of  the  dependent  variable  for  a  unit 
change  of  the  particular  independent  variable  while  the  remaining 
independent  variables  are  held  statistically  constant.  It  is, 
however,  possible  to  standardize  the  partial  regression 
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TABLE  1   :  MULTIPLE  REGRESSION  EQUATIONS 
MINIMUM  OZONE  LEVEL  PREDICTION 
(MAY  1988  to  SEPTEMBER  1989) 


The  Regression  Equation  is  of  the  form: 
Y  =  Bo  +  Bi  *Xi  +  B2  *Xa  +  Ba  *X3  +  B*  *X« 


Y   -  Maximum  Ozone  Level  (PPM) 
Xi   -  Maximum  Temperature   (C) 
Xa   -  Average  Wind  Speed   (Km/hr) 
Xa   -  Average  Cloud  Opacity  (  10*^  "  ■  ) 
X4   -  Previous  Day's  Maximum  Ozone  (PPM) 
Recorded  at  Simcoe  (22071) 


TORONTO 

(31104) 

[LRT] 


TORONTO 
(31104) 


SCARBOROUGH 

(33033) 

[LRT] 


SCARBOROUGH 
(33033) 


NORTH  YORK 
(34020) 
[LRT] 


NORTH  YORK 
(34020) 


N 

Bo 
Bi 
B} 
Ba 

B4 

R 

R» 

SE 

Y 

Xi 

Xi 

Xa 

X4 


35 

-71.798 

5.166 

-1.026 

-0.524 

0.152 

0.805 

0.648 

15.410 

72.818 

29.230 

14.153 

1.923 

72.659 


41 

■54.829 

3.931 

-0.894 

-1.625 

0.115 

0.679 

0.461 

19.777 

55.133 

26.380 

13.128 

1.895 

55.967 


36 

-59.189 

4.985 

-1.330 

-0.679 

0.125 

0.698 

0.487 

17.515 

76.111 

29.229 

13.952 

1.880 

72.756 


42 

■23.973 

3.333 

-1.088 

-2.525 

0.099 

0.653 

0.427 

19.556 

56.000 

26.350 

12.886 

1.862 

56.563 


30 

-66.524 

5.335 

-1.047 

-1.902 

0.109 

0.685 

0.470 

21.621 

77.268 

29.461 

14.231 

2.177 

74.146 


35 

-27.924 

3.381 

-0.949 

-1.413 

0.051 

0.696 

0.485 

12,941 

57.241 

26.938 

15.257 

1.801 

57.931 


N  -  Number  of  Observations 

Bo  -  Multiple  Regression  Coefficient  (intercept) 

Bi  -  Multiple  Regression  Coefficient  (for  Xi  ) 

Ba  -  Multiple  Regression  Coefficient  (for  Xa  ) 

Ba  -  Multiple  Regression  Coefficient  (for  Xa  ) 

B4  -  Multiple  Regression  Coefficient  (for  X*  ) 

R  -  Correlation  Coefficient 

R»  -  R-Squared 

SE  -  Standard  Error  of  the  Estimate  (PPB) 

Y  -  Mean  Measured  Ozone  (PPB) 

Xi  -  Mean  Maximum  Temperature  (C) 

Xa  -  Mean  Average  Wind  Speed  (Km/hr) 

Xa  -  Mean  Average  Cloud  Opacity 

X4  -  Mean  Previous  Day's 

Measured  Ozone  at  Simcoe  (PPB) 
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coefficients  with  new  values  represented  as  bj  's,  derived  as 
follows  (King,  1969): 

bj  =  Bj    S  (2) 

(To 

where  Oj  and  0^  are  the  standard  deviations  of  the  independent 
and  dependent  variables  respectively.   Since  they  are 
dimensionless  items,  the  bj  's  can  be  compared  directly  as 
measures  of  the  relative  importance  of  the  different  variables  in 
accounting  for  the  variation  in  the  dependent  variable.   The  bj  's 
for  each  site  is  found  in  Table  2. 

In  all  cases,  the  most  important  variable  is  air 
temperature.   For  the  LRT  cases,  air  temperature  is  about  3  times 
more  important  in  the  regression  equation  than  the  next  most 
important  variable,  wind  speed.   The  importance  of  wind  speed  in 
the  regression  for  non-LRT  cases,  however,  is  not  much  less  than 
that  of  air  temperature.   Cloud  opacity  and  the  previous  day's 
ozone  play  the  smallest  role  in  the  regression  equation,  however, 
their  inclusion  provides  some  improvement  in  the  prediction. 

The  overall  significance  of  the  multiple  regression  can  be 
tested  by  an  analysis  of  variance.   This  test  is  equivalent  to 
testing  the  null  hypothesis  that  in  the  population  from  which  the 
observations  have  been  sampled,  there  is  no  significant 
regression  of  the  dependent  on  the  independent  variables.   The 
test  consists  of  determining  the  F  ratio  for  each  case,  that  is 
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TABLE  2  :   STANDARDIZED  REGRESSION  COEFFICIENTS 
AND  ANALYSIS  OF  VARIANCE  TEST 
(MAY  1988  to  SEPTEMBER  1989) 


TORONTO 

(31104) 

[LRT] 


TORONTO 
(31104) 


SCARBOROUGH 

(33033) 

[LRT] 


SCARBOROUGH 
(33033) 


NORTH  YORK 

(34020) 

[LRT] 


NORTH  YORK 
(34020) 


1 

0.699 

0.587 

0.603 

0.567 

0.650 

0.505 

<1 

-0.197 

-0.347 

-0.242 

-0.488 

-0.189 

-0.370 

13 

-0.037 

-0.156 

-0.045 

-0.282 

-0.120 

-0.133 

14 

0.144 

0.139 

0.111 

0.040 

0.098 

0.062 

"l 

3.633 

3.136 

3.592 

3.069 

3.661 

3.151 

12 

5.389 

8.163 

5.410 

8.093 

5.422 

8.211 

"3 

1.981 

2.013 

1.997 

2.019 

1.903 

2.020 

'* 

26.516 

25.428 

26.687 

24.764 

27.011 

25.622 

îr 

28.048 

20.983 

29.687 

18.032 

30.015 

21.103 

14.  375* 

6.279* 

8.653- 

6.088* 

10.161* 

6.622* 

bi  =  Standardized  Regression  coefficient 

bi  =  Standardized  Regression  coefficient 

ba  =  Standardized  Regression  coefficient 

b4  =  Standardized  Regression  coefficient 

cfi  =  Standard  Deviation  for  Xi 

Oi  =  Standard  Deviation  for  Xa 

03  =  Standard  Deviation  for  X3 

64  =  Standard  Deviation  for  X4 

<5t  =  Standard  Deviation  for  Y 


(X.  ) 
(X.  ) 
(X3) 

(X4  ) 


F  =  F  -  Ratio 


Significcuit  at  5  %  level  of  confidence 


the  ratio  of  the  mean  square  of  the  sxaa   of  squares  of  the 
regression  over  the  mean  square  of  the  sum  of  squares  of  the 
residual.   If  the  computed  F  exceeds  the  tabled  value  for  a 
selected  level  of  confidence,  then  the  null  hypothesis  is 
rejected. 

The  F  ratio  is  found  at  the  bottom  of  Table  2 .   In  all 
cases,  the  null  hypothesis  can  be  rejected  at  5%  level  of 
confidence. 

The  effect  of  LRT  from  the  industrial  United  States  can  be 
determined,  at  any  particular  site,  from  the  difference  in  the 
predicted  ozone  levels  using  identical  values  for  the  independent 
variables  in  both  the  LRT  and  non-LRT  regression  equations.   This 
difference  was  determined  by  using  the  above  equations  and  all 
combinations  of:  air  temperature  being  equal  to  23,  27,  30  and  34 
degrees  C,  wind  equal  to  0,  5,  10  and  20  km/hr.,  cloud  cover 
equal  to  0,  2,  and  4  tenths,  and  the  previous  day's  ozone  equal 
to  40,  60,  80,  100  and  120  ppb.   The  results  showed  a  minimiim 
increase  due  to  LRT  to  be  24.2%  (13.7  ppb),  a  maximum  increase  of 
48.7%  (38.1  ppb)  and  an  average  of  31.7%  (25.9  ppb).   However,  it 
can  be  expected  that  the  effect  of  LRT  to  increase  ozone 
concentrations  would  be  even  greater  than  indicated  in  this 
analysis  which  does  not  take  into  account  the  difference  in  the 
previous  day's  ozone  levels  during  LRT  and  non-LRT  days.   The 
mean  of  the  previous  day's  ozone  is  15  ppb  larger  during  LRT  than 
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non-LRT  days  (see  Table  1)  and  all  of  the  regression  equations 
are  positively  correlated  with  the  previous  day's  ozone. 

5.    MODEL  TEST 

The  multiple  regression  equation,  as  a  predictor  of  peak 
level  ozone,  was  tested  by  using  independent  data  sets  from  May 
to  September  1987  for  the  Toronto  and  Scarborough  sites. 
Measured  peak  ozone  levels  were  compared  against  predicted  peak 
ozone  levels  by  using  the  appropriate  regression  equation,  the 
relevant  meteorology  and  the  measured  previous  day's  ozone  at 
Simcoe.   Figures  1  to  4  show  graphs  of  measured  versus  predicted 
ozone  levels  for  each  of  LRT  and  non-LRT  cases  at  Toronto  and 
Scarborough  while  Figure  5  shows  the  total  measured  versus 
predicted  graph. 

Willmott  (1982),  citing  several  examples  showing  that  real 
significance  and  statistical  significance  do  not  exhibit  a  close 
correspondence,  has  proposed  measures  of  'systematic'  and 
'unsystematic'  error  to  examine  the  types  of  differences  and  an 
'index  of  agreement'  to  investigate  the  relative  size  of 
differences.   Because  a  model  should  explain  most  of  the  major 
trends  present  in  the  observed  record,  it  is  important  to  know 
what  proportion  of  RMSE  is  systematic  (RMSEa  )  and  what  proportion 
is  unsystematic  (RMSEu ) .   In  this  fashion,  with  respect  to  a 
'good'  model,  (RMSE.  )=  should  approach  zero  while  (RMSEu  )=  should 
approach  RMSE.   The  'index  of  agreement'  (D)  was  proposed  to 
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yield  a  measure  of  the  relative  size  of  differences;  this  index 
is  intended  to  be  a  descriptive  measure  which  is  both  relative 
and  bounded  (0  <.  D  <.  1)  and  can  be  widely  applied  in  order  to 
make  cross-comparisons  between  models .   As  the  agreement  between 
the  observed  and  predicted  variables  increases,  D  approaches 
1.0  .   Table  3  shows  the  equations  used  to  calculate  model 
performance . 


Table  4  contains  the  model  performance  statistics .   It 
should  first  be  noted  that  the  calculated  and  measured  population 
means  and  standard  deviations  are  very  similar,  in  fact  in  no 
case  is  the  difference  between  the  two  means  greater  than  2  ppb. 

The  MAE  and  RMSE  indicate  a  mean  difference  between  the 
predicated  and  observed  values  of  approximately  13  ppb.   With  the 
unsystematic  portion  of  the  RMSE  being  around  10  ppb,  this 
indicates  that  most  of  the  RMSE  is  unsystematic  and  it  need  not 
be  assvimed  that  there  is  a  missing  independent  variable  in  the 
regression  analysis. 

The  last  column  of  Table  4  contains  the  model  performance 
statistics  for  the  combined  "computed  versus,  measured"  matrix. 
This  is  included  to  increase  the  size  of  the  population  to  be 
'large'  (N>30)  since  the  general  utility  of  such  indices  as  RMSE 
is  questionable  as  N  becomes  small  (Willmott,  1982).   The  high 
value  of  the  Willmott  Index  of  Agreement  (D)  (0.853)  indicates  a 
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TABLE   3  :    Statistical  Measures   of  Model   Performance 
(after  Willmot,    19  82) 


I.      SUMMARY  STATISTICS: 


1.     Mean  Observed; 


2.      Mean  Predicced; 


3.      Scandard  Deviation: 
(Observed) 


-1  ^ 
0     =     N        Z     0. 

1=1      ' 

N 
Z 

i=l 

N 


p    =    n"-^  Z    P 


n"^  Z   (0.  -  0)^ 
L         i=l      ^ 


4.      Scandard  Deviation: 
(Predicted) 


-1  ^      -2 
N   Z  (P.  -  P)^ 

L    i=i   ^ 


5.   Intercept  (a)  and  Slope 
(b)  of  the  Least  Squares 
Regression: 


P.      =     a  +  b  0. 

1  i 


II.      DIFFERENCE  STATISTICS: 


6.     Mean  Biased  Error: 


MBE     =     N   ^   Z    (P.    -  0.) 
i=l      "  ^ 


7.      Mean  Absolute  Error: 


MAE  =  N  ^  Z  If.  -  0. 
i=I   ^    " 


8.   Root  Mean  Square  Error:       RMSE  = 


-1  ^         2 
N   Z  (P.  -  0.) 

1=1 


9.   Root  Mean  Square 
Systenatic  Error: 


RMSE, 


N  ^  Z  (P.  -  0.)^ 

1=1 


10.   Root  Mean  Square 

Unsystematic  Error: 


RMS 


^ 


-1  ^       -  2 
N  '  Z  (P.  -  P.) 

i=l  ^    ^  . 


11.   Index  of  Agreement: 


D  =   1  - 


I  (?.  -  0  )y  z  (|p  -  ol+lo  -  ol)2 

i=l         /i=l    "■        ^ 


.  0  <  D  <  1 
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TABLE  4   :  OZONE  LEVEL  PREDICTION 

MODEL  PERFORMANCE  STATISTICS 


TORONTO        TORONTO       SCARBOROUGH    SCARBOROUGH     TOTAL 
[LRT]  [LRT] 


N 

12 

16 

11 

15 

54 

0 

69.000 

48.312 

73.273 

55.866 

60.092 

P 

69.112 

50.291 

74.004 

56.517 

61.033 

s. 

19.821 

17.647 

22.865 

17.410 

21.198 

Sp 

17.102 

13.424 

16.714 

10.907 

16.915 

a 

0.603 

0.528 

0.470 

0.422 

0.604 

b 

27.512 

24.757 

39.547 

32.953 

24.755 

MBE 

0.112 

1.978 

0.731 

0.550 

0.940 

MAE 

11.872 

10.557 

15.256 

10.556 

11.806 

RMSE 

13.928 

12.498 

16.816 

12.479 

13.794 

RMSE. 

7.537 

8.295 

11.572 

9.747 

8.376 

RMSE. 

11.713 

9.348 

12.201 

7.791 

10.959 

R» 

0.488 

0.483 

0.413 

0.453 

0.573 

D 

0.821 

0.816 

0.764 

0.763 

0.853 

N  -  Niimber  of  Observations 

0  -  Average  Observed  Ozone 

P  -  Average  Predicted  Ozone 

So  -  Observed  Standard  Deviation 

Sp  -  Predicted  Standard  Deviation 

a  -  Slope  of  Regression  Line  (OBS  vs  PRED) 

b  -  Intercept  of  Regression  Line 

MBE  -  Mean  Bias  Error 

MAE  -  Mean  Absolute  Error 

RMSE  -  Root  Meem  Square  Error 

RMSE.  -  Systematic  Portion  of  RMSE 

RMSE.  -  Unsystematic  Portion  of  RMSE 

R'  -  R-Squared 

D  -  Wilmott  Index  of  Agreement 
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reasonably  good  agreement  between  the  predicted  and  measured  peak 
ozone  levels . 

6.    8 -HOUR  AVERAGE  OZONE  LEVELS 

The  final  part  of  this  study  deals  with  the  relationship 
between  the  peak  hourly  ozone  level  and  the  8-hour  average  on  the 
same  day.   The  data  set  used  for  the  determination  of  the 
regression  equations  is  used  here  to  determine  the  linear 
regression  between  the  8 -hour  average  (Y)  and  the  peak  hourly 
average  (Xo )  by  the  equation: 

Y  =  Bo  +  Bi  *Xi  (2) 

Table  5  shows  the  linear  regression  equation  for  each  site 
and  for  each  of  LRT  and  non-LRT  days.   From  the  table  we  can  see 
a  high  level  of  correlation  with  an  average  correlation 
coefficient  (R)  of  0.95  and  an  average  coefficient  of 
determination  (R^  )  of  0.89  .   Using  these  linear  equations  it  is 
possible  to  estimate  the  8 -hour  average  associated  with  a  peak 
hourly  ozone  level  of  80  ppb  (Moderate  air  quality)  and  120  ppb 
(Poor  air  quality) .   These  values  are  listed  in  Table  5  and  are 
equal  to,  on  average,  62  and  91  ppb,  respectively.   Figures  6  to 
11  show  graphs  of  the  8-hour  average  ozone  level  versus  the 
maximum  ozone  level  for  each  of  LRT  and  non-LRT  cases  at  Toronto, 
Scarborough  and  North  York  while  Figure  12  shows  the  total  8-hour 
average  versus  maximum  ozone  graph. 
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TABLE  5   :  LINEAR  REGRESSION  EQUATIONS 

8  HOUR  AVERAGE  OZONE  LEVEL  PREDICTION 
(MAY  1988  to  SEPTEMBER  1989) 

The  Regression  Equation  is  of  the  form: 
Y  =  Bo  +  Bi  *Xi 


Y   -  8  Hour  Average  Ozone  (PPM) 
Xi   -  Peak  Hourly  Ozone     (PPM) 

TORONTO    TORONTO    SCARBOROUGH    SCARBOROUGH    NORTH  YORK    NORTH  YORK 
(31104)    (31104)      (33033)        (33033)       (34020)       (34020) 
[LRT]  [LRT]  [LRT] 


N 

Bo 

Bi 

R 

R' 

SE 

Y 

Xi 


35 
0.495 
0.772 
0.985 
0.918 
5.993 
56.152 
72.818 


MOD.   62.225 

LEVEL 

POOR   93.135 

LEVEL 


41 

10.177 

0.631 

0.943 

0.889 

5.430 

41.379 

55.133 

60.657 

85.897 


36 
1.886 
0.766 
0.943 
0.889 
6.296 
59.494 
76.111 

63.166 

93.806 


42 
4.479 
0.682 
0.964 
0.930 
5.220 
41.941 
56.000 

59.039 

86.319 


30 

35 

-4.196 

5.210 

0.855 

0.714 

0.945 

0.948 

0.893 

0.899 

8.512 

5.244 

62.383 

45.610 

77.268 

57.241 

64.204 

62.330 

98.404 

90.810 

N 

Bo   - 

Bi 

R 

R» 

SE   - 

Y 

Xi 
MOD.  - 
LEVEL 
POOR  - 
LEVEL 


Number  of  Observations 

Linear  Regression  Coefficient  (intercept) 

linear  Regression  Coefficient  (for  Xi  ) 

Correlation  Coefficient 

R-Squared 

Standard  Error  of  the  Estimate  (PPB) 

Mean  8  Hour  Average  Ozone  (PPB) 

Mean  Peak  Hourly  Ozone  (C) 

Value  of  8  Hour  Average  Calculated 

for  a  Peak  Ozone  Level  of  80  PPB 

Value  of  8  Hour  Average  Calculated 

for  a  Peak  Ozone  Level  of  120  PPB 
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7.    SUMMARY  AND  CONCLUSIONS 

This  study  has  determined  that  the  levels  of  ozone  in 
Metropolitan  Toronto  increases  on  LRT  days .   These  increases  in 
ozone  averaged  approximately  32%  (26  ppb)  and  peaks  as  high  as 
49%  (38  ppb)  were  noted.   These  results  are  compared  with  Wolff 
and  Lioy  (1978),  who  concluded  by  a  similar  means  that  the 
contribution  of  external  sources  to  a  Connecticut  location  during 
high  ozone  days,  was  as  much  as  7  5  ppb,  and  Prior  et  al.  (1981) 
reported  the  LRT  contribution  to  St.  Louis  was  36  ppb. 

Multiple  regression  equations  were  determined  relating  the 
peak  hour  ozone  levels  with  the  same  day  maximxim  temperature,  the 
average  daytime  windspeed,  the  cloud  opacity  at  Pearson 
International  Airport  and  the  peak  hour  ozone  concentration  at 
the  Simcoe  rural  station  for  the  previous  day,  for  LRT  and  non- 
LRT  days  in  1988  and  1989.   Data  for  1987  were  used  to  validate 
the  statistically  determined  relationships .   These  equations  have 
proven  to  be  reasonably  good  predictors  of  peak  ozone  levels  in 
the  Metropolitan  Toronto  area. 

The  study  also  indicated  that  Ontario's  80  ppb  one-hour 
average  regulatory  limit  cannot  be  used  as  a  surrogate  for  an 
eight-hour  average  health  standard.   The  World  Health 
Organization  Guidelines  recommend  50  to  60  ppb  for  an  eight-hour 
average  level  of  ozone  for  the  protection  of  human  health. 
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Similar  studies  could  be  carried  out  in  the  future  to 
determine  regression  equations  for  peak  hourly  ozone  for  other 
cities  of  Southern  Ontario,  for  both  LRT  and  non-LRT  days.   The 
regression  ecpiations  can  be  used  to  predict  the  peak  ozone  levels 
and  the  associated  AQI  for  each  location.   A  good  suggestion  by  a 
Ministry  Scientist  is  to  use  the  previous  day's  maximiim  ozone 
level  at  the  C.N.  Tower  elevated  station  as  a  prediction 
variable.   This  could  be  tested  in  a  future  study. 
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APPENDIX  A 


Maximum  Observed  Ozone 
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APPENDIX  B 
Number  of  Ozone  Potential  Days 
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